Light Quark Resonances in pp Annihilations at 5.2 GeV/c 
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Data from the Fermilab E835 experiment have been used to study the reaction pp — > 777771-° at 5.2 
GeV/c. A sample of 22 million six photons events has been analyzed to construct the Dalitz plot 
containing ~ 80fc 777777° events. A partial wave analysis of the data has been done. Six /j-states 
decaying into 7777 and five aj-states decaying into 7771-° are identified in the mass region ~ 1.3 and 
2.4 GeV, and their masses, widths and spins are determined by maximum likelihood analysis of 
the data. Two fo states are identified with the popular candidates for the lightest scalar glueball, 
/o(1500) and /o(1710). 



PACS numbers: 14.40.Cs 13.75.Cs 13.25.Jx 
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I. INTRODUCTION 

Antiproton-proton annihilations leading to final states 
containing light-quark structures constitute very attrac- 
tive and powerful tools for the study of the spectroscopy 
of qq mesons, glueballs, and hybrids. Indeed, major con- 
tributions to light-quark spectroscopy have been made 
by the study of pp interactions with stopped and low en- 
ergy antiprotons Q]. In this paper we report a study of 
the isoscalar and isovector light quark resonances popu- 
lated in the reaction pp — > 777771"° in the collision of 5.2 
GeV/c antiprotons with protons at rest. The reaction 
is very selective. It can only populate qq states with 



jPC 



■,++ 



isoscalar, or /-states decaying into 7777, 



and isovector, or a-states decaying into 7777. The 7777 
states may contain admixtures of low lying ++ and 2 ++ 
glueballs, and the 7777° states may include states with 
non-qq J PC , such as the 1 h hybrids. 

In the spectroscopy of light-quark mesons there remain 
many open questions. Theoretical calculations 2] for 
states with angular momentum L < 4 predict twelve 
/j-states, and six aj-states in the mass range 1 GeV- 
2.5 GeV. In addition, two glueball states, with 1 = 0, 
J PC = ++ ,2 ++ , and at least one manifestly exotic qqg 
hybrid with 1 = 1, J PC = 1 h are predicted in this mass 
range. In different experiments many fo and f% states 
have been claimed, but PDG'04 Q considers only four 
fo states, four fa states, and five 00,2,4 to be established 
well enough to be included in the meson summary list. 
In the present investigation we present new determina- 
tions of masses and J PC of several of the established fj 
and a j states, and provide evidence for several previously 
unconfirmed states. 



II. DATA AND EVENT SELECTION 

In the present analysis data from Fermilab pp annihi- 
lation experiment E835 is used. The stored and cooled 



antiprotons circulating in the Fermilab Antiproton Accu- 
mulator intersect a hydrogen cluster jet gas target, and 
the reaction products are detected in a detector system 
which surrounds the interaction region. Full details of 
the detector and its performance can be found elsewhere 
4]. Photons are detected in the lead glass central detec- 
tor in the angular range 10.6° < < 70.0° and 27r in <f), 



with an energy resolution of 



<te{E) 
E 



+ 0.04. 



^E(GeV)) 

Luminosity was measured by detecting recoil protons at 
~ 87.5°. The data were collected at p momenta of 5.2 
GeV/c corresponding to center of mass energy range of 
Ecms = 3.409 - 3.418 GeV with a total luminosity of 
UApb- 1 . 

The reaction studied was 



pp — > Tr riri, 77" 



n 



27, 77 -> 27 
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Twenty two million six photon events were recorded in 
the full detector, of which 3.25 million were fully con- 
tained in the fiducial volume, had no photons in the for- 
ward calorimeter, and met the energy- momentum (4C) 
constraint with probability > 10%. Of these, 139,500 
events passed 777777° seven constraint fit with CL > 10%. 
These events were then subjected to three 'anticuts'. 
Events which met the hypotheses for the 'contaminant' 
final states (7r 7r 7r°), (7777° tt ) or (777777) with CL > 2% 
were removed, and the final sample of 83,400 events was 
used for the Dalitz plot analysis. The Dalitz plot and its 
7777 and 7777° projections are shown in Fig. ^ These mass 
distributions differ completely from the phase space dis- 
tributions which are predicted to be essentially flat over 
the entire 1.0 — 2.6 GeV/c 2 mass region for both M{jirf) 
and M(?77r ). Visible structures appear to correspond in 
the 7777 diagonal to / (1500), /o(1710), and a complex of 
states in the mass region of ~ 2000 - 2250 MeV. The 
visible horizontal and vertical 7777° bands appear to cor- 
respond to a (980) and a state around 1300 MeV. It is 
important to keep in mind that the 'anticuts' described 
above remove events and introduce distortions in specific 
regions of the Dalitz plot. The effect of the 37r° 'anti- 
cut' was large, but mostly confined to the diagonal edge 
of the Dalitz plot in the /2(1270) region. The effect of 
the 7r 7r°77 'anticut' was smaller and also largely confined 
to the / 2 (1270) and / (1370) regions. The effect of 377 
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'anticut' was small, and was mainly in the vicinity of 
02(1320). 

Monte Carlo simulation were done to estimate feed- 
down contributions from 7 and 8 photon reaction with 
one and two missing photons, respectively, which remain 
in our final event selection. It was found that the contri- 
butions of ljtjt: (5.4 ± 0.2) %, and w7T°7r° (5.2 ± 0.2)%, 
were mostly confined to the diagonal edge of the Dalitz 
plot, in the region of / 2 (1270) and / (1370). The contri- 
butions of 47T° was estimated as (2.3 ± 0.1)%, and was 
mostly confined to the region around a 2 (1320). The over- 
all efficiency of our i]i]ir° selection was estimated by MC 
simulation to be (4.26 ± 0.1)%. Taking account of the 
feed-down contributions in the total yield, we obtain the 
total cross section corresponding to the 777771"° events in 
the Dalitz plot, a^n ) = 150 ± l(stat) ± 2{syst) nb. 



III. FORMALISM FOR DALITZ-PLOT 
ANALYSIS 

The reaction pp — > 777777° is described in the isobar 
model as a two step process, pp — > fj + tt°, or aj + 77, 
followed by fj — > 7777, and aj —> rjir . In principle both 
the production and the decay of fj and aj should be in- 
dependently taken into account. In practice, this is not 
possible. At an incident antiproton momentum of 5.2 
GeV/c, pp annihilation involves angular momentum up 
to ~ 7 in the initial state. This makes it essentially im- 
possible to use the full helicity formalism for the analysis 
of our data. Instead, we use the decay formalism develo- 
ped and successfully used in several Dalitz plot analyses 
of the Crystal Barrel annihilation data at 1.94 GeV/c 
ULEi III- I n addition, we take account of the production 
amplitudes in an average way by parameterizing them as 
a function of even powers of cos 9, where 9 is the reso- 
nance production angle with respect to the p direction, 
p(9) = 1 + Y^i=i Pi c °s 2z 0. The final amplitude for the 
formation and decay of a resonance with spin J, projec- 
tion A = J x , mass mo and width Tq is the given by 



A) =p(9)G x e i ^F J (q)- 



(a,/3) 



im T r , 



(1) 



where G\ and the 5\ are the magnitude and the phase 
of the complex coupling constant, Yj(a,{3) are spherical 
harmonics which are functions of angles a and /3, the po- 
lar and azimuthal angles of the decay in the rest frame 
of the resonance. The denominator is that for the rela- 
tivistic Breit Wigner with energy dependent width, 



r, 



r f^M) I ( V'foo) 

V m J V m o 



(2) 



Here Fj(q) are the Blatt-Weisskopf barrier factors for the 
resonance break-up momentum q (subscripts denote 
values at the resonance mass, mo), evaluated for a radius 
of 1 fermi. The combined intensities of two interfering 



resonances 1, and 2 in the Dalitz plot are given by 

wi,2= ExUxa\ 2 + \A X , 2 \ 2 + 2c x ^(A XA Al 2 )], (3) 

where the interference coefficient c\ ranges from (no co- 
herence) to ±1 (maximum coherence). We wish to note 
here that when two resonances interfere the difference be- 
tween their phase angles, A^^t) also becomes important. 
In considering interferences a certain pragmatic compro- 
mise has to be made. With 11 different resonances in our 
final result (see Table|TT| , if all possible two state interfer- 
ences, including self-interference were taken into account, 
one would need (60 x A multiple) values of c\. This is 
impossible and also unnecessary. In the present work a 
large number of possible interferences were studied, but 
only those were retained in the final analysis which im- 
proved the log-likelihood by more than 20 per additional 
degree of freedom. The negative log-likelihood of a MC 
fit to the data was defined as 



InL = Nln( 



M 



N 



)-<E h 



.DATA 



), 



(4) 



where wf ATA refer to the N events in the data, and wf 10 
refer to the M events in the MC fit. 

With this definition an increase of A(lnL) = (0.5 x r) 
with increase of r degrees of freedom corresponds to a one 
standard deviation (A% 2 /d.o.f. = 1) improvement in the 
fit. For making the final fits mo, To and J are manually 
varied, p{9) are separately fit, and the free parameters for 
the fit are G\, S\ and c\. As detailed below, in practice 
much more demanding criteria for A(lnL) were used. 

The ingredients of fitting the Dalitz plot are masses, 
widths, J PC , the production amplitude parameters 
p(9), and resonance-resonance interferences including 
self-interference. We take these into account succes- 
sively. We start with a basic set of six resonances. 
These include those which are identifiable as clear en- 
hancements in the Dalitz plot, /o(1500), /o(1710), and 
/2(2150) as representative of the broad enhancements in 
the M(rm) ~ 2000 - 2200 MeV region, and a 2 (1320). 
In this basic set we also include a 2 (1700) and £14 (2040) 
which are not apparent as identifiable enhancements, but 
which are considered well established and whose decay in 
7771° has been observed 0]. 

The initial fit with this set of resonances was done at 
the simplest level. Their masses and widths were fixed 
at PDG values neither the the angular distribution 
of the production amplitudes, nor resonance-resonance 
interference was included, and only the ten G\ were fit- 
ted by maximizing —InL. The fit was quite poor, with 
lnL/#par. = -9,520/10. 

As a first step in improving the fit, it was decided to 
take account of the fact that at 5.2 GeV/c the production 
amplitudes are far from isotropic. It was found that their 
angular distributions could be satisfactorily parameter- 
ized as a third order polynomial in cos 2 (9) . This added 18 
parameters in the fit which improved by AlnL/ Aftpar. = 
5,261/18. Inclusion of the strong self-interference of 
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FIG. 1: (a) Dalitz plot of rjr/Tv . (b) and (c) show rjri and nir mass projections. The positions of the resonances used in the 
initial fit are indicated at their PDG mass values in MeV. The histograms show the relatively poor fit obtained with this initial 
set of resonances, as described in the text. 




FIG. 2: (a) The best fit for the Dalitz Plot, and the corresponding invariant mass projections, (b) r/r/, and (c) r/ir . Data are 
shown with points. Histograms show the best fit using the final set of resonances. 



Resonance 


AlnL 1 A#param. 




Interference 


AlnL/ A#param. 


CO 


A(5 ) 


Cl 


A(5r) 


/o(2020) 


1464/1 




a 4 (2040) x a 4 (2040) 


363/2 


-0.51 





-1.00 





/o(2100) 


158/1 




a 4 (2040) x a 2 (1320) 


168/4 


0.83 


0.06 


1.00 


0.41 


/ 2 (2340) 


878/2 




a 4 (2040) x / 2 (2150) 


155/4 


-0.48 


1.29 


-1.00 


1.36 


ao(1450) 


297/1 




a 4 (2240) x a 4 (2240) 


254/2 


-0.80 





-0.83 





a 4 (2240) 


267/2 




a 2 (1320) x /o(2100) 


161/2 


-1.00 


0.17 







TABLE I: Resonances and interferences added in the second iteration. c\ and A(<5a) (A = 0, 1) refer to the interference 
coefficients and the phase differences (in radians) between the interfering resonances. 



a 4 (2040) led to an additional AlnL/A#par. = 363/2, 
bringing the net —InL/^par. = 15,144/30. The resul- 
ting fit is shown in Fig. \i\(b) and (c), for the projec- 
tions of invariant masses M(rjT]) and M(i]ir ), respec- 
tively. Both are rather poorly fit, and it is obvious that 
additional resonances and interferences must be consid- 
ered. 

The additional resonances input include /q(2020), 
/ (2100), / 2 (2340), and a (1450), which are not con- 
sidered by PDG as firmly established, and a new state 
a4(2240). As shown in Tabic [I] the log-likelihood im- 



provements due to the inclusion of these resonances is in 
all cases more than 130 per added parameter. Further, 
only those interferences are included which increase log- 
likelihood by more than ~ 40 per added parameter. Fig. 
121 shows the much improved fit to the mass projections. 

For the finally chosen resonances two dimensional op- 
timization of the mass and width was done for different 
J PC assignments. The final result of the optimization 
was -lnL/#par. = 19894/44. 
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FIG. 3: Results for A| — lnL\ optimization of mass and width: (a) contour plot for /o(1710), (b) A| — lnL\ variation for mass 
and width for /,/(1710), J PC = 0++ and 2++, (c) same for /j(1500) 
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FIG. 4: Results for A| - lnL\ optimization for (a) /o(2020), (b) /o(2100), and / 2 (2150). 
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FIG. 5: Results for Aj - lnL\ optimization for (a) / 2 (2340), (b) aj(1320), and a (1450). 



IV. RESULTS AND DISCUSSION 

Fig. [21 shows the overall results of the best fit to 
the Dalitz plot and the mass projections for M(jjrj) and 
M{r)Tr). Our final results are summarized in Table II, 
along with the corresponding results from PDG 0. We 
wish to note that the relative intensity fractions listed 
in Table II are related to, but not the same as relative 



branching fractions. This is partly due to the approx- 
imations inherent in the formalism used in our partial 
wave analysis, and partly because the anticuts we have 
imposed in event selection have different effects in dif- 
ferent parts of the Dalitz plot. Also, interferences have 
been omitted in evaluating these fractions, which there- 
fore do not add to 100%. The last column in the table 
lists the change in log-likelihood if a particular resonance 
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is dropped from the fit, and the parameters of all other 
resonances are reoptimized. This provides information 
about the relative importance of a resonance which is 
complementary to that given by its relative intensity. We 
now briefly discuss results shown in Figs. 3 — 6 for indi- 
vidual resonances. 

Our procedure for obtaining the best fit values of the 
parameters of a resonance is illustrated in Fig. EJa,b) 
with the results for /j(1710). The | — lnL\ is calculated 
for specific values of mass, width, and J PC of /o(1710). 
The results are fitted with smooth curves to determine 
the best fit values of M and T corresponding to the max- 
ima in | - lnL\. Fig. GJb) shows that the J PC = 0++ 
solution is clearly preferred. The corresponding contour 
plot for M/r is shown in Fig. Efa). The same proce- 
dure was followed for /j(1500), whose results are shown 
in Fig. OUc). Again J PC = ++ is clearly preferred by 
A(— InL) = 46. The results for other resonances are 
shown are shown in Figs. 01 03 and[fj] In all but one case 
the | — InL | differences between the different spin choices 
are large, and the spin assignment with the largest |— lnL\ 
(shown with solid points) is preferred; the others, (shown 
with open points) are rejected. For a j (2040) (Fig. 6), 
although the | — lnL\ values for J PC = 2 ++ and 4 ++ arc 
comparable, no maxima are found for either M or T for 
J PC = 2 ++ , and it is also rejected. 



A. The Isoscalar States fo and /2. 

As already mentioned, Fig. [3] clearly shows that both 
/o(1500) and / (1710) have J PC = 0++. The J PC = 
++ assignment is preferred over 2 ++ by AlnL = 58 for 
/o(1500), and AlnL = 46 for / (1710). As listed in Table 
ITT1 for /o(1500) our mass is ~ 30 MeV smaller than the 
PDG average, but the width is identical. For /q(1710) 
our mass is ~ 30 MeV larger, and the width is also ~ 30% 
larger. We note that the Crystal Barrel has reported that 
while they do observe /o(1500), they do not find any ev- 
idence for /o(1710) excitation in the pp — > 777777° reaction 
in pp annihilation at either 0.9 GeV/c 8] or 1.94 GeV/c 



On the other hand, WA102 |jj reports excitation of 
both /o(1500) and /o(1710) in central production with 
decay into r\r\ . We find that in our pp — > 777777° measure- 
ment at 5.2 GeV/c, in the Dalitz plot the ratio of the 
intensities (/o(1710) -► W)/(/o(1500) -> r)r)) ~ 2.2. We 
note that Amsler et al. |8| report this ratio to be < 0.25 
(90% CL) for production in pp annihilation at 0.9 GeV/c. 
Hopefully, our result provides a new input in the contin- 
uing discussion about the scalar glueball admixtures in 
these two states. 

Although the observation of the isoscalars Jo (2020) 
and /o(2100) has been reported before, PDG'04 consi- 
ders them unconfirmed. We find / (2020) and / (2100) 
with masses and widths listed in the Table [H] to be 
strongly populated with nearly equal intensity. Once 
again the ++ assignments are preferred for both over 
2 ++ assignments by AlnL = 23 and 22, respectively. 
In a global analysis of the Crystal Barrel data for mo- 
menta 900-1900 MeV/c Anisovich et al. [uj identify 
these two scalars decaying into 7777 as having mass/ width 
of (2005±30)/(305±50) MeV and (2105±15)/(200±25) 
MeV, which are in almost exact agreement with ours. 
We also note that a combined enhancement of these two 
states was reported as X(2100) decaying into 7777 by the 
earlier Fermilab E760 experiment [ll|. Its J PC was not 
determined. 

Although ten /2-states have been reported in our mass 
region by different experiments, the PDG regards only 
four, /£(1525), / 2 (2010), / 2 (2300) and / 2 (2340) as suffi- 
ciently well established to be included in the meson sum- 
mary list. None of these except /2(1525) have been ob- 
served to decay into 7777. We find no evidence for the 
excitation of /2(1525), although it is possible that its 
resolution is rendered impossible due to the presence of 
nearby stronger /o(1500). The InL plot in Fig. [^c) sug- 
gests that / 2 (1525) population is very small, if not zero. 
Similarly, we find no evidence for / 2 (2010). We confirm 
the existence of / 2 (2150) which has been omitted by PDG 
from its summary list, but which was previously observed 
to decay into 7777 by WA102 Our mass and width are 
in good agreement with both WA102 and the PDG ave- 
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Mass(MeV) 




Width(MeV) 


Rel. Int. 


If removed 


Resonance 


PDG'04 


Present 


J PC 


PDG'04 


Present 


% 


AlnL 
A#par. 


/o(1500) 


1507 ± 5 


1473 ± 5 


0++ 


109 ±7 


108 ±9 


2.4 


207/1 


/o(1710) 


1714 ±5 


1747 ± 5 


0++ 


140 ± 10 


188 ±13 


5.2 


521/3 


/o(2020)* t 


1992 ± 16 


2037 ± 8 


0++ 


442 ± 60 


296 ± 17 


10.9 


295/3 


/o(2100)* t 


2103 ± 7 


2105 ± 8 


0++ 


206 ± 15 


236 ± 14 


7.9 


249/3 


/ 2 (2150)* 


2156 ± 11 


2170 ± 6 


2++ 


167 ± 30 


182 ± 11 


6.5 


209/6 


/ 2 (2340)* t 


2339 ± 55 


2350 ± 7 


2++ 


31 Q+ 81 


218 ±16 


4.8 


124/2 


a 2 (1320) 


1318 ± 1 


1327 ± 2 


2++ 


111 ± 2 


128 ±4 


6.5 


2335/14 


a (1450) 


1474 ± 19 


1477 ± 10 


0++ 


265 ± 13 


267 ±11 


3.9 


418/5 


a 2 (1700)* 


1732 ± 16 


1702 ± 7 


2++ 


194 ± 40 


417 ±19 


20.6 


1185/4 


a 4 (2040) 


2010 ± 12 


2004 ± 6 


4 ++ 


353 ± 40 


401 ± 16 


20.1 


738/13 


a 4 (2237)* 




2237 ± 5 


4 ++ 




291 ± 12 


20.7 


823/4 



TABLE II: Results for masses, widths and J of light quark resonances decaying into rjrj and r/n as seen in pp — > 7r°7777 at 
5.2 GeV/c. All errors are statistical only and correspond to A| — lnL\ = 0.5 or A% 2 = ±1<j. Asterisks* mark states presently 
omitted from meson summary list of PDG'04. Daggers^ mark f-states which have not been seen in r\r] before. For explanation 
of the last two columns, see text. 



rage. We also observe clear evidence for /2(2340). Etkin 
et al. 01 had claimed this as a doublet, /2(2300) and 
/2(2340), in its <fxj) decay. We do not find any evidence 
for such a close mass doublet. 

We note that we do not find any evidence for the tensor 
glueball candidate £(2230), or / 2 (2230) decaying into rjij. 
This is consistent with the result of the Crystal Barrel 
search for the same [la ■ 



B. The Isovector States ao, a 2 and a 4 . 

Our fit to the r/ir mass spectrum (see fig. |2Jc)) is 
admittedly poorer compared to the isoscalar state in fig 
HJb) . It also presents a certain vexing problem concerning 
the strong enhancement observed at mass ~ 1320 MeV, 
which is followed by a dip in the vicinity of 1450 MeV. 
The a2(1320) is a well established tensor resonance with 
PDG mass/width, M/T = (1318.3 ± 0.6)/(105 - 111) 
MeV. It is strongly populated in many reactions, and 
its decays into several channels, including rjir , are well 
documented 0. A scalar a (1450) with M/T = (1470 ± 
25)/ (265 ± 30) MeV decaying into i]tt° was reported by 
CB ^2 in. pp annihilation at rest, and also in K^K^ |15| . 

On the other hand, OBELIX [lfjj has reported a scalar 
with M/T = (1290 ± 10)/(80± 5) in pp annihilation de- 
caying into K^Kg, along with a nearly factor five weaker 
excitation of a 2 (1320) than what follows from the PDG 
averages. We have examined the possible existence of 
a (1290) of OBELIX. As shown in Fig. GJb) and Ta- 
bic [n] wc obtain a good fit to our data with a 2 (1320) 
and ao(1450), with the best fit M/T parameters for 
both, a 2 : M/T = (1327 ± 2)/(128± 4) MeV and a : 
M/T = (1477 ± 10)/(267± 11), both in excellent agree- 
ment with the PDG'04 values. Paradoxically, an even 



better fit (not shown) is obtained with do (1290) with 
M/T = (1314 ± 2)/(133 ± 4) MeV and a (1450) with 
M/T = (1404±9)/(222±18) MeV. With a 66% (or 8.3a) 
larger width for ao(1290) than that reported by OBELIX 
(133 ± 4 MeV versus 80 ± 5 MeV) it is unlikely that this 
solution represents their resonance. However, this still 
poses the question whether both 02 (1320) and ao(1290) 
can coexist. We have studied this question at length, 
considered additional interferences, and conclude that if 
both resonances exists, it is not possible to resolve them. 
We list the results for the J — 2 assignment in Table 
[n] because of the weight of numerous final states, inclu- 
ding 7771", in which a 2 (1320) has been observed, and not 
ao(1290). Nevertheless, we must consider the ambiguity 
between a 2 (1320) and ao(1290) as unresolved. 

We confirm 04 (2040) with parameters in excellent 
agreement with those in PDG'04. We confirm a 2 (1700) 
but find it to have much larger width than PDG'04. We 
find a strong resonance, 04 (2237) decaying into 7771, with 
M/T = (2237 ± 5)/(291± 15) MeV. A J PC = 4++ reso- 
nance decaying into 7777 and tj'tt was earlier reported by 
Anisovich et al. |l7|. 

Additional resonances which were tried in the fit 
were a (980), / (1370), #(1525), / 2 (1640), / 2 (1910), 
7Ti(1400), tti(1600), / (1710 - 1790), / 2 (1870 - 1910). 
None of these obtain significant contributions (> 0.5%) 
to the intensity. As mentioned earlier, the resonances 
near the boundaries of the Dalitz Plot, particularly 
ao(980) and /o(1370), suffer from the effect of the 'an- 
ticuts'. No convergence was obtained when the reported 
hybrids ?Ti(1400) and tti(1600) were included in the fits. 

The authors wish to thank the E835 Collaboration for 
allowing us the use of their data. This work was sup- 
ported by the US Department of Energy. 
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